VLT/SPHERE robust astrometry of the HR8799 planets at milliarcsecond-level accuracy Orbital architecture analysis with PyAstrOFit

HR8799 is orbited by at least four giant planets, making it a prime target for the recently commissioned Spectro-Polarimetric High-contrast Exoplanet REsearch (VLT/SPHERE). As such, it was observed on five consecutive nights during the SPHERE science verification in December 2014. We aim to take full advantage of the SPHERE capabilities to derive accurate astrometric measurements based on H-band images acquired with the Infra-Red Dual-band Imaging and Spectroscopy (IRDIS) subsystem, and to explore the ultimate astrometric performance of SPHERE in this observing mode. We also aim to present a detailed analysis of the orbital parameters for the four planets. We report the astrometric positions for epoch 2014.93 with an accuracy down to 2.0 mas, mainly limited by the astrometric calibration of IRDIS. For each planet, we derive the posterior probability density functions for the six Keplerian elements and identify sets of highly probable orbits. For planet d, there is clear evidence for nonzero eccentricity ($e \simeq 0.35$), without completely excluding solutions with smaller eccentricities. The three other planets are consistent with circular orbits, although their probability distributions spread beyond $e = 0.2$, and show a peak at $e \simeq 0.1$ for planet e. The four planets have consistent inclinations of about $30\deg$ with respect to the sky plane, but the confidence intervals for the longitude of ascending node are disjoint for planets b and c, and we find tentative evidence for non-coplanarity between planets b and c at the $2 \sigma$ level.Comment: 23 pages, 14 figure

Spectral cube extraction for the VLT/SPHERE IFS: Open-source pipeline with full forward modeling and improved sensitivity

We present a new open-source data-reduction pipeline to reconstruct spectral data cubes from raw SPHERE integral-field spectrograph (IFS) data. The pipeline is written in Python and based on the pipeline that was developed for the CHARIS IFS. It introduces several improvements to SPHERE data analysis that ultimately produce significant improvements in postprocessing sensitivity. We first used new data to measure SPHERE lenslet point spread functions (PSFs) at the four laser calibration wavelengths. These lenslet PSFs enabled us to forward-model SPHERE data, to extract spectra using a least-squares fit, and to remove spectral crosstalk using the measured lenslet PSFs. Our approach also reduces the number of required interpolations, both spectral and spatial, and can preserve the original hexagonal lenslet geometry in the SPHERE IFS. In the case of least-squares extraction, no interpolation of the data is performed. We demonstrate this new pipeline on the directly imaged exoplanet 51 Eri b and on observations of the hot white dwarf companion to HD 2133. The extracted spectrum of HD 2133B matches theoretical models, demonstrating spectrophotometric calibration that is good to a few percent. Postprocessing on two 51 Eri b data sets demonstrates a median improvement in sensitivity of 80% and 30% for the 2015 and 2017 data, respectively, compared to the use of cubes reconstructed by the SPHERE Data Center. The largest improvements are seen for poorer observing conditions. The new SPHERE pipeline takes less than three minutes to produce a data cube on a modern laptop, making it practical to reprocess all SPHERE IFS data.Comment: 17 pages, 11 figures. Software available at: https://github.com/PrincetonUniversity/charis-de

Searching for faint companions with VLTI/PIONIER. II. 92 main sequence stars from the Exozodi survey

The Exozodi survey aims to determine the occurrence rate of bright exozodiacal discs around nearby main sequence stars using infrared interferometry. Although the Exozodi survey targets have been carefully selected to avoid the presence of binary stars, the results of this survey can still be biased by the presence of unidentified stellar companions. Using the PIONIER data set collected within the Exozodi survey, we aim to search for the signature of point-like companions around the Exozodi target stars. We use both the closure phases and squared visibilities collected by PIONIER to search for companions within the ~100 mas interferometric field of view. The presence of a companion is assessed by computing the goodness of fit to the data for a series of binary models with various separations and contrasts. Five stellar companions are resolved for the first time around five A-type stars: HD 4150, HD 16555, HD 29388, HD 202730, and HD 224392 (although the companion to HD 16555 was independently resolved by speckle interferometry while we were carrying out the survey). In the most likely case of main sequence companions, their spectral types range from A5V to K4V. Three of these stars were already suspected to be binaries from Hipparcos astrometric measurements, although no information was available on the companions themselves so far. In addition to debiasing the statistics of the Exozodi survey, these results can also be used to revise the fraction of visual binaries among A-type stars, suggesting that an extra ~13% A-type stars are visual binaries in addition to the ones detected in previous direct imaging surveys. We estimate that about half the population of nearby A-type stars could be resolved as visual binaries using a combination of state-of-the-art interferometry and single-aperture imaging, and we suggest that a significant fraction of these binaries remains undetected to date.Comment: Accepted for publication in A&

Discretized aperture mapping with a micro-lenses array for interferometric direct imaging

Discretized Aperture Mapping (DAM) appears as an original filtering technique easy to play with existing adaptive optics (AO) systems. In its essential DAM operates as an optical passive filter removing part of the phase residuals in the wavefront without introducing any difficult-to-align component in the Fourier conjugate of the entrance pupil plane. DAM reveals as a new interferometric technique combined with spatial filtering allowing direct imaging over a narrow field of view (FOV). In fact, the entrance pupil of a single telescope is divided into many sub-pupils so that the residual phase in each sub-pupil is filtered up to the DAM cut-off frequency. DAM enables to smooth the small scale wavefront defects which correspond to high spatial frequencies in the pupil plane and to low angular frequencies in the image plane. Close to the AO Nyquist frequency, such pupil plane spatial frequencies are not well measured by the wavefront sensor (WFS) due to aliasing. Once bigger than the AO Nyquist frequency, they are no more measured by the WFS due to the fitting limit responsible for the narrow AO FOV. The corresponding image plane angular frequencies are not transmitted by DAM and are useless to image small FOVs, as stated by interferometry. That is why AO and DAM are complementary assuming that the DAM cut-off frequency is equal to the AO Nyquist frequency. Here we describe the imaging capabilities when DAM is placed downstream an AO system, over a convenient pupil which precedes the scientific detector. We show firstly that the imaging properties are preserved on a narrow FOV allowing direct imaging throughout interferometry. Then we show how the residual pupil plane spatial frequencies bigger than the AO Nyquist one are filtered out, as well as the residual halo in the image is dimmed

A near-infrared interferometric survey of debris-disc stars. IV. An unbiased sample of 92 southern stars observed in H-band with VLTI/PIONIER

Context. Detecting and characterizing circumstellar dust is a way to study the architecture and evolution of planetary systems. Cold dust in debris disks only traces the outer regions. Warm and hot exozodiacal dust needs to be studied in order to trace regions close to the habitable zone. Aims. We aim to determine the prevalence and to constrain the properties of hot exozodiacal dust around nearby main-sequence stars. Methods. We search a magnitude limited (H < 5) sample of 92 stars for bright exozodiacal dust using our VLTI visitor instrument PIONIER in the H-band. We derive statistics of the detection rate with respect to parameters such as the stellar spectral type and age or the presence of a debris disk in the outer regions of the systems. We derive more robust statistics by combining our sample with the results from our CHARA/FLUOR survey in the K-band. In addition, our spectrally dispersed data allows us to put constraints on the emission mechanism and the dust properties in the detected systems. Results. We find an over-all detection rate of bright exozodiacal dust in the H-band of 11% (9 out of 85 targets) and three tentative detections. The detection rate decreases from early type to late type stars and increases with the age of the host star. We do not confirm the tentative correlation between the presence of cold and hot dust found in our earlier analysis of the FLUOR sample alone. Our spectrally dispersed data suggest that either the dust is extremely hot or the emission is dominated by the scattered light in most cases. The implications of our results for the target selection of future terrestrial planet finding missions using direct imaging are discussed.Comment: 20 pages, 16 figures, 4 tables, updated references and minor changes to the text, accepted for publication in A&

Fundamental limitations of high contrast imaging set by small sample statistics

In this paper, we review the impact of small sample statistics on detection thresholds and corresponding confidence levels (CLs) in high contrast imaging at small angles. When looking close to the star, the number of resolution elements decreases rapidly towards small angles. This reduction of the number of degrees of freedom dramatically affects CLs and false alarm probabilities. Naively using the same ideal hypothesis and methods as for larger separations, which are well understood and commonly assume Gaussian noise, can yield up to one order of magnitude error in contrast estimations at fixed CL. The statistical penalty exponentially increases towards very small inner working angles. Even at 5-10 resolution elements from the star, false alarm probabilities can be significantly higher than expected. Here we present a rigorous statistical analysis which ensures robustness of the CL, but also imposes a substantial limitation on corresponding achievable detection limits (thus contrast) at small angles. This unavoidable fundamental statistical effect has a significant impact on current coronagraphic and future high contrast imagers. Finally, the paper concludes with practical recommendations to account for small number statistics when computing the sensitivity to companions at small angles and when exploiting the results of direct imaging planet surveys.Comment: 12 pages, 10 figures, accepted to Ap

Deep imaging of β Pictoris at L’: asymmetries in the disc and constraints on planets

We demonstrate here the potential of VLT/NaCo at L’ (3.8μm) in Angular Differential Imaging (ADI). We detect the β Pictoris disc above 5σ between 0.4’’ and 3.8’’, after combining data from 7 epochs or 3.1 years. To avoid the smearing of the planet due to its orbital motion within this time span (green arrow), we subtracted the planet from the 7 datasets. We used a star subtraction technique based on PCA (Soummer et al. 2012) and corrected for ADI biases by iterating. We implemented a forward modeling approach to constrain the dust distribution

Discovery of a Companion Candidate in the HD169142 Transition Disk and the Possibility of Multiple Planet Formation

We present L' and J-band high-contrast observations of HD169142, obtained with the VLT/NACO AGPM vector vortex coronagraph and the Gemini Planet Imager, respectively. A source located at 0".156+/-0".032 north of the host star (PA=7.4+/-11.3 degrees) appears in the final reduced L' image. At the distance of the star (~145 pc), this angular separation corresponds to a physical separation of 22.7+/-4.7 AU, locating the source within the recently resolved inner cavity of the transition disk. The source has a brightness of L'=12.2+/-0.5 mag, whereas it is not detected in the J band (J>13.8 mag). If its L' brightness arose solely from the photosphere of a companion and given the J-L' color constraints, it would correspond to a 28-32 MJupiter object at the age of the star, according to the COND models. Ongoing accretion activity of the star suggests, however, that gas is left in the inner disk cavity from which the companion could also be accreting. In this case the object could be lower in mass and its luminosity enhanced by the accretion process and by a circumplanetary disk. A lower mass object is more consistent with the observed cavity width. Finally, the observations enable us to place an upper limit on the L'-band flux of a second companion candidate orbiting in the disk annular gap at ~50 AU, as suggested by millimeter observations. If the second companion is also confirmed, HD169142 might be forming a planetary system, with at least two companions opening gaps and possibly interacting with each other.Comment: Accepted to ApJL, see also Biller et al. 201

SPHERE IRDIS and IFS astrometric strategy and calibration

We present the current results of the astrometric characterization of the VLT planet finder SPHERE over 2 years of on-sky operations. We first describe the criteria for the selection of the astrometric fields used for calibrating the science data: binaries, multiple systems, and stellar clusters. The analysis includes measurements of the pixel scale and the position angle with respect to the North for both near-infrared subsystems, the camera IRDIS and the integral field spectrometer IFS, as well as the distortion for the IRDIS camera. The IRDIS distortion is shown to be dominated by an anamorphism of 0.60+/-0.02% between the horizontal and vertical directions of the detector, i.e. 6 mas at 1". The anamorphism is produced by the cylindrical mirrors in the common path structure hence common to all three SPHERE science subsystems (IRDIS, IFS, and ZIMPOL), except for the relative orientation of their field of view. The current estimates of the pixel scale and North angle for IRDIS are 12.255+/-0.009 milliarcseconds/pixel for H2 coronagraphic images and -1.75+/-0.08 deg. Analyses of the IFS data indicate a pixel scale of 7.46+/-0.02 milliarcseconds/pixel and a North angle of -102.18+/-0.13 deg. We finally discuss plans for providing astrometric calibration to the SPHERE users outside the instrument consortium.Comment: 12 pages, 6 figures, 3 table

Apocenter pile-up and arcs: a narrow dust ring around HD 129590

Observations of debris disks have significantly improved over the past decades, both in terms of sensitivity and spatial resolution. At near-infrared wavelengths, new observing strategies and post-processing algorithms allow us to drastically improve the final images, revealing faint structures in the disks. These structures inform us about the properties and spatial distribution of the small dust particles. We present new $H$-band observations of the disk around HD 129590, which display an intriguing arc-like structure in total intensity but not in polarimetry, and propose an explanation for the origin of this arc. Assuming geometric parameters for the birth ring of planetesimals, our model provides the positions of millions of particles of different sizes to compute scattered light images. We demonstrate that if the grain size distribution is truncated or strongly peaks at a size larger than the radiation pressure blow-out size we are able to produce an arc quite similar to the observed one. If the birth ring is radially narrow, given that particles of a given size have similar eccentricities, they will have their apocenters at the same distance from the star. Since this is where the particles will spend most of their time, this results in a "apocenter pile-up" that can look like a ring. Due to more efficient forward scattering this arc only appears in total intensity observations and remains undetected in polarimetric data. This scenario requires sharp variations either in the grain size distribution or for the scattering efficiencies $Q_\mathrm{sca}$. Alternative possibilities such as a wavy size distribution and a size-dependent phase function are interesting candidates to strengthen the apocenter pile-up. We also discuss why such arcs are not commonly detected in other systems, which can mainly be explained by the fact that most parent belts are usually broad.Comment: Accepted for publication in A&A, abstract shortene